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We present the next-to-next-to-leading order (aaj) corrections to the first moment of the polarized virtual 
photon structure function g2{x,Q^, P^) in the kinematical region Q^ ^ P^ ^ A^ in QCD. We find that the 
aa^ corrections are about 3% of the sum of the leading (a) and the next-to- leading (aas) contributions, when 
g2 = 30 ~ lOOGeV^ and P^ = 3GeV2. 



1. Introduction 

The investigation of the photon structure has 
been an active field of research both theoretically 
and experimentally |1I2I3I4| . Also there has been 
growing interest in the study of the spin struc- 
ture of photon. In particular, the first moment of 
the polarized photon structure function gj has at- 
tracted attention in connection with its relevance 
to the QED and QCD axial anomaly |5I6I7I8I9| . 
The polarized photon structure functions can be 
measured from two-photon processes in the po- 
larized e~^e~ collider experiments as shown in 
Fig. 1, where —Q^{—P^) is the mass squared of 
the probe (target) photon. 

For a real photon (P^ = 0) target, there ex- 
ists only one spin-dependent structure function 
g1{x,Q'^). The QCD analysis for gJ was per- 
formed in the leading order (LO) ^D] and in 
the next-to-leading order (NLO) |11I12| . In the 
case of a virtual photon target (P^ ^ 0) there 
appear two spin-dependent structure functions, 
gl{x,Q^,P'^) and g'^ {x , Q"^ , P'^) . The former has 
been investigated up to the NLO in QCD by the 
present authors in |13ll4| . and also in the second 
paper of ^^l- In Ref.|13| the structure function 
gj{x, Q^,P^) was analysed in the kinematical re- 



Figure 1. Deep inelastic scattering on a virtual 
photon in e+ e^ collision. 



gion 

A^ < p2 < : 



(1) 
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where A is the QCD scale parameter. The ad- 
vantage to study the virtual photon target in this 
kinematical region is that we can calculate struc- 
ture functions by the perturbative method with- 
out any experimental data input ^B], which is 
contrasted with the case of the real photon tar- 
get where in the NLO there exist nonperturbative 
pieces. 

In this talk we focus on the photon structure 
function gJ and, in particular, on the sum rule 
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and present the result of the NNLO ((aa^)) cor- 
rections to the first moment of gJ{x,Q^,P^) in 
the case of the kinematical region Q • 

2. The sum rule of gj 

The polarized structure function gJ of the real 
photon satisfies a remarkable sum rule 



dx gJ{x,Q^) =0 . 



(2) 



In fact, the authors of Ref. [H| showed that the 
sum rule (0), holds to all orders in perturbation 
theory in both QED and QCD. 

When the target photon becomes off-shell, i.e., 
P^ 7^ 0, the first moment of the corresponding 
photon structure function gJ {x , Q'^ , P'^) does not 
vanish any more. Indeed, for the case A^ <^ P^ -^ 
Q^, the first moment has been calculated up to 
the NLO (0(aa^)) in QCD as follows |7in| : 



J dxgUx,Q\P^) 
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with /3o = 11 — 2n//3 being the one-loop QCD 
/3 function. Here a ~ e^/Air and as{Q^) — 
'g^{Q^)/4:TT are the QED coupling constant and 
the QCD running coupling constant, respectively, 
and Ci is the electromagnetic charge of the active 
quark with flavour i in the unit of proton charge 
and 71/ is the number of active quarks. Note that 
the r.h.s. of (j^J does not involve any experimen- 
tal data input. This is because the nonpertur- 
bative pieces, which take part in the real photon 
target case, can be neglected in the kinematical 
region |^ and the whole twist-two contributions 
to gj{x, Q^, P^) can be computed by the pertur- 
bative method. It is also noted that the first term 
in the square brackets of the r.h.s. of © results 
from the QED triangle anomaly while the second 
term comes from the QCD triangle anomaly. 



3. The NNLO (aal) corrections 

Now we calculate the NNLO (ctaf) corrections 
to the first moment of gjlx, Q^,P^). First recall 
that for the operator product expansion (OPE) 
of two electromagnetic (and thus gauge- invariant) 
currents, only gauge-invariant operators need be 
included with their renormalization basis [Ifij . 
Since there is no gauge-invariant twist-two gluon 
and photon operators with spin one, we need to 
consider only quark operators, i.e., the flavour 
singlet Rg and nonsinglet i?^5 axial currents, as 
follows: 

Rs = ?^7"75 IV^, Rns = V^7"75(Qch-(e')l)^, (4) 
where 1 is an n/ x n/ unit matrix and Q^f^ is the 
square of the Uf x Uf quark-charge matrix. 

Then the first moment of gJ{x,Q^,P^) is ex- 
pressed as 



dxgJix,Q^,P') 



(5) 



l2 /.,2 



,a)(J{p)\Rsi^^^)hip)) 

Here Cs and Cns are the coefficient functions 
corresponding to the currents Rg and R%[g, re- 
spectively, and (7(p)|i?,;(/i^)|7(p)) with i = S, NS 
are the photon matrix elements of these quark ax- 
ial currents. We choose the renormalization point 
H at //2 =p2_ 

The coefficient functions are given by 
a{QyP',g{P'),a) 



exp 



9(P' 



SiQ^ 



dg 



Pig') 



a{l,g{Q^),a), (6) 



where i — S,NS and 7^(5) is the anomalous di- 
mension of the quark axial current Rf and (3{g) is 
the QCD /3- function. We expand 74(5) in powers 
of g as 



7^(5) = 7 



(0) g 



9 



167r 



IGtt' 



:f 



2+^^^( 



(7) 



Since the nonsinglet quark axial current R%g is 
conserved in the massless limit, it undergoes no 
renormalization, and thus we have 



(0) _ (1) _ (2, _ 



= ••• = 



(8) 
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On the other hand, the singlet axial current Rg 
has a non-zero anomalous dimension 75(5) due to 
the axial anomaly. At the one-loop level, we know 
7^°) = 0. The two-loop [181 and three-loop P| 
results are: 



7i^^ = 



UCpUf, 
■CfCa 



(2M 



V 3 



ZQCl 



V 



° n 2 



(9) 
(10) 



.(2) 



with Cf = I and Ca — 3. The result 7^ ' 
obtained in the MS scheme with a definition of 
the 75 -matrix as 

7/^75 = ^'^f^parYl"!^, (11) 

being used. 

The P function has been calculated up to the 
four-loop level in the MS scheme 20 21 . For nu- 
merical analysis, we will use later the QCD run- 
ning coupling constant as{Q^) where the results 
up to the three- loop level are taken care of j22| . 
But for the present we only need the /? function 
up to the two-loop level: 

dg 



M 



dfi 



P{9) 



-/3o 



9 



/3i 



9 



167r2 ''^'(167r2)2 
with the SUci'i) value 
2 



+ 



(12) 



/3o 

Pi 



11 



V 



38 
102- -nf 



(13) 



Using Eqs.Q and ((T^ we obtain 



exp 



aiP' 



aiQ^ 



dg 



ns{9') 

P{9') 



I 

8/3, 



7s ' (as{Q'^) 
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We already have results of the singlet coeffi- 
cient function Cs(l,5(Q^), a) calculated up to 
the two-loop level |23l24j and the nonsinglet coef- 
ficient function CNs(.^ig{Q'^)^ ot) up to the three- 
loop level (23 • Both were calculated in the MS 
scheme and with the definition of the 75-matrix 
in idJ. Taking up to the two- loop level, we have 



Cs(l,g(Q'),«)/(e2) 
r2l 



1 



3^ a.(Q' 



-Cf 



32 



Cf 



23^ ,1 13, 

Y^Ca + (2C3 + ^)n/^ 



«.(g' 



C^s(i,g(0' 



3^ as{Q^) 



(15) 
(16) 



+Cf{^^Cf 



.32 



23 
16 



Ca 



1 



v)( 



a.(Q') 



where C3 is the Riemann zeta-function (^3 = 

1.202056903 •••)• 

For — p2 = p2 -;^ ^2^ ^]-^g photon matrix el- 
ements of the quark axial currents can be calcu- 
lated perturbatively and they are expressed in the 
form as 



{^{p)\R,{^Ji^^p^)\^{p)) = -A. 



An 



with i 



A, 



4°)-H 



5, /V5, and 



167r2 



A, 



(1) 



167r2 



A 



(2)_ 



(17) 
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The leading terms Ai ' are connected with Adler- 
Bell-Jackiw anomaly and are given bv|26| 



A 



(0) 



= -12n/(e') , 



,(0) 



A;^>, = -I2nf{{e^) ~ {e'Y) . 



(19) 
(20) 



In our talk given at RADCOR05, we reported 
that for the next-leading and next-next-leading 



terms A] ' and A] \ 



we have A 



A 



(2) 



A 



(2) 
NS 



1) 
S 



A 



(1) 

NS 



due to the nonrenormalization 



theorem |17j for the triangle anomaly. It is true 
that we obtain 



4(1) _ 4(1) _ /i(2) _n 



(21) 



(2) 



but we found after the workshop that Ag has 
a non- vanishing value. For the calculation of 
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Ag , we need to evaluate the three- loop Feynman 
graphs. Instead we resort to the Adler-Bardeen 
theorem for the axial current Jl^{— Rg in Q), 



d.J^- ' ""' 



^•^5 



167r2 2 






(22) 



which holds in all orders in as, and we evaluate, 
in the MS scheme, the matrix element 

eA''""(0|T[A^a,A?^^(-p)A,(p)]|0}amputated , (23) 

in the two-loop level, where A^ and A^ are giuon 
and photon fields, respectively. Details will be 
reported elsewhere I^T). We found 



Af^2An,{e-)C,'^C4- 



2 V 3 



8C; 



(24) 



M 



Putting all these equations into Eq.jSJ with 
P^, we finally obtain the NNLO [aal) cor- 



rections to the first moment of g1{x, Q^, P^): 
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where the third to 8th terms are the NNLO con- 
tributions. In the case of nj = 4, for an example, 
we have 



f dxgl{x,Q\P^) 
Jo 



3a as(Q^) ocsiP'^) 

= <^ 0.4198 - 0.1235 ^^ ' - 0.2963^^^ — '- 

-0.02731 ( ^^ M +0.01185 ^^ ' ^^ — - 



TT / TT TT 

?2U 2I 



-1.153 



(^)- 



(26) 



To estimate the sizes of the NLO (aag) and 
NNLO (aa^) corrections compared to the LO (a) 
term, we take, for instance, Q^ — 30 and 100 
GeV^, and P^ = 3 GeV^. The corresponding 
values of as are obtained from Ref.{22^. We get 
as{Q^ = 30GeV2) = 0.2048, a^Q^ = lOOGeV^) = 
0.1762, and as{P^ = 3GeV^) = 0.3211. The re- 
sults are given in Tabled 



4. Summary 

We have investigated the next-to-next-to- 
leading order (aa^) corrections to the first mo- 
ment of the polarized virtual photon structure 
function gJ{x,Q^,P^) in the kinematical region 
Q^ > P^ > A^ in QCD. AU the necessary infor- 
mation on the coefficient functions and anoma- 
lous dimensions corresponding to the quark ax- 
ial currents has been already known, except for 
the three-loop-level photon matrix element (the 
finite term) of the fiavour singlet quark axial cur- 
rent Rg . Instead of evaluating the relevant three- 
loop Feynman diagrams, we resort to the Adler- 
Bardeen theorem for the axial current H22(l . Then 
calculation reduces to the one in the two-loop 
level. We evaluate in effect the two-loop diagrams 
for the photon matrix element of the gluon oper- 
ator (ESI- 

The aa^ corrections are found to be about 3% 
of the sum of the leading (a) and the next-to- 
leading {aas) contributions, when Q^ = 30 '^ 
lOOGeV^ andp2 = 3GeV^ 



NNLO corrections to the first moment of gj{x, Q"^, P^) 

Table 1 

The NLO and NNLO contributions relative to LO in the first moment of gj{x, Q^,P^) 







QVGeV" 


P' 


-/GeV" 


LO 


NLO 


NNLO 


NNLO/(LO+NLO) 


Uf-- 


= 3 


30 




3 




-0.0816 


-0.0250 


-0.0272 






100 




3 




-0.0766 


-0.0234 


-0.0254 


nf~- 


= 4 


30 




3 




-0.0913 


-0.0288 


-0.0317 






100 




3 




-0.0886 


-0.0287 


-0.0315 


"^f-- 


= 5 


30 




3 




-0.0986 


-0.0309 


-0.0343 






100 




3 




-0.0977 


-0.0319 


-0.0354 
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